Fidelity decay in trapped Bose-Einstein condensates 
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The quantum coherence of a Bose-Einstein condensate is studied using the concept of quantum 
fidelity (Loschmidt echo). The condensate is confined in an elongated anharmonic trap and subjected 
to a small random potential such as that created by a laser speckle. Numerical experiments show 
that the quantum fidelity stays constant until a critical time, after which it drops abruptly over a 
single trap oscillation period. The critical time depends logarithmically on the number of condensed 
atoms and on the perturbation amplitude. This behavior may be observable by measuring the 
interference fringes of two condensates evolving in slightly different potentials. 

PACS numbers: 03.75.Gg, 05.45.Mt 



Introduction. — Ultracold atom gases have been at the 
center of intensive investigations since the first realiza- 
tion of an atomic Bose-Einstein condensate (BEC) in the 
mid-1990s. Applications of BECs include the possibil- 
ity of revisiting standard problems of condensed-matter 
physics [l| , by making use of periodic optical lattices that 
mimic the ionic lattice in solid-state systems. More re- 
cently, BECs have been used to study quantum transport 
in disordered systems (another long-standing problem in 
condensed matter physics), by using a laser speckle to 
create a disordered potential. For instance, it was shown 
that the free expansion of a BEC is restrained, or even 
completely suppressed, in the presence of a disordered 
potential UHf (see also Q for a theoretical analysis), an 
effect akin to Anderson localization in solids [f|. 

These problems are often approached by studying the 
interference pattern of two or more BECs that are re- 
leased from the optical trap, expand freely, and even- 
tually interact with each other [aHH 9]. Experiments 
show high-contrast matter-wave interference fringes, thus 
revealing the coherent nature of Bose-Einstein conden- 
sates. Surprisingly, recent experiments have shown high- 
contrast fringes even for well separated BECs, whose 
phases are totally uncorrelated [1 01 ] . 

Because the interference pattern depends on the phases 
of the condensates, the fringes should be sensitive to per- 
turbations that strongly affect the phase, but weakly af- 
fect the motion. Therefore, when two condensates are 
subjected to a random potential (such as that generated 
by a laser speckle [TJ EI) before interfering, we expect 
a reduction in the contrast of the interference fringes, 
which should depend both on the amplitude of the ran- 
dom potential and on the time during which it has been 
in contact with the condensates. The purpose of this Let- 
ter is to propose a theoretical procedure to quantify this 
loss of coherence and to suggest a possible experimental 
realization. 

In order to estimate the coherence and stability of a 
quantum system [TU, one can compare the evolution of 
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the same initial condition in two slightly different Hamil- 
tonians, Hi — Hq + SHi and Hi = H0 + 6H2, where Hq is 
the unperturbed Hamiltonian, and 8H± t 2 are small per- 
turbations characterized by the same amplitude, same 
wavelength spectrum, but different phases. The quan- 
tum fidelity at time t is then defined as the square of the 
scalar product of the wavefunctions evolving with H\ and 
H 2 respectively: F(t) — {{i^Hx (*)] 4>h 2 (*))P- This proce- 
dure is sometimes referred to as the 'Loschmidt echo', as 
it is equivalent to evolving the system forward in time 
with H±, then backward with Hi, and using the fidelity 
to check the accuracy of the time-reversal. 

Virtually all theoretical investigations of the 
Loschmidt echo consider one-particle systems evolving in 
a given (usually chaotic) Hamiltonian. Several regimes 
have been described in the past. For perturbations that 
are classically weak but quantum-mechanically strong, 
the fidelity decay is exponential, with a rate independent 
on the perturbation and given by the classical Lya- 
punov exponent of the unperturbed system (l4|. This 
behavior has been confirmed by numerical simulations 

, ■ For weaker perturbations, the decay rate is still 
exponential, but perturbation-dependent (Fermi golden 
rule regime). For still weaker perturbations, the decay 
is Gaussian (perturbative regime) [HI]. For integrable 
systems, other types of decay (notably algebraic) have 
been observed [17] . Finally, a perturbation-independent 
regime, though with Gaussian decay, was also observed 
in experiments [Lsl ] . 

In a previous work [l^| , we have applied the concept of 
Loschmidt echo to a system of many electrons interact- 
ing through their self-consistent electric field. The nu- 
merical results showed that the quantum fidelity remains 
equal to unity until a critical time, then drops suddenly 
to much lower values. A similar result was also obtained 
for a classical system of colliding hard spheres [2C| . This 
effect is probably related to the nonlinearity introduced 
by the interactions between particles. Therefore, BECs 
should constitute an ideal arena to determine whether 
such behavior is typical of many-body quantum systems. 

Model. — The dynamics of a BEC is accurately de- 
scribed, in the mean-field approximation, by the Gross- 
Pitaevskii equation (GPE). We considered a cigar-shaped 
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condensate, where the transverse frequency of the con- 
fining potential is much larger then the longitudinal fre- 
quency, to± 3> lo z . In this case, a one-dimensional (ID) 
approximation can be used, and the GPE reads as: 



h 2 d 2 ijj 
2m d z 2 



+ V(z)T/j + g 1D N A \^\ 2 ij; = H ^j . (1) 



m^ + Kz^/L 2 ^), 
harmonic oscillator 



Here, \ijj\ 2 dz = 1, Na is the number of condensed 
atoms, giD — 2ahio±_ is the ID effective coupling con- 
stant, and a is the 3D scattering length. The confining 
potential contains a small quartic component, which can 
be realized optically [21(, V(z) = ~™ Ji :j 
where L\ w = (h/mcjz) 1 / 2 is the 
length. 

We choose the parameters of the experiment described 
in Ref. [2], where N A = 10 5 atoms of 87 Rb (a = 5.7 nm) 
arc confined in a quasi-lD trap with lo z /2tt — 24.7 Hz and 
uj±/2tt — 293 Hz. In the simulations, we normalize time 
to Ug, space to L\ lo = 2.16 /im, and energies to Hu> z . 
The dimensionless ID coupling constant is then: gio = 
9iD/(L\ 10 huj z ) = 0.063. The quartic coefficient, which 
will be crucial to excite a sufficiently complex nonlinear 
dynamics, is taken to be K = 0.05. With the above 
parameters, the half-length of the condensate is roughly 
24 /im. 

The random perturbation SH can be realized in prac- 
tice using a laser speckle [lH fl2| . In our simula- 
tions, we model the random potential through the sum 
of a large number of uncorrelated waves: 8H/Huj z — 
e Ylj=N ■ cos(2ttz/ Xj + ctj), where e is the amplitude 
of the perturbation, Aj's are the wavelengths, and ctj's 
are random phases. The smallest wavelength present in 
the random potential is A m j n = 2Li 10 = 4.32 /im, which 
is consistent with the experimental correlation length, 
a z = 5 fim. [2j . The wavelength spectrum of the pertur- 
bation (i.e. the values of N m i n and iV max ) affects only 
weakly the behavior of the fidelity: therefore, we will fo- 
cus our analysis on the dependence of the fidelity on the 
amplitude e. 

In order to compute the quantum fidelity, we proceed 
as follows: (i) first, we prepare the condensate in its 
ground state without perturbation; (ii) then, we sud- 
denly displace the anharmonic trap by a distance Az 
(a few micrometers); (iii) finally, we solve numerically 
the time-dependent GPE ((T|) with the perturbed Hamil- 
tonian Ho + 5H . Step (iii) is performed for N = 11 
uncorrelated realizations of the random potential, thus 
yielding N evolutions of the wavefunction, ipj(t). We 
then use all possible combinations to compute the par- 
tial fidelities Etj (t) = \(i>i{t)\i } j(t))f- There are of course 
AI = Nl/(N — 2)!2! = 55 independent combinations, 
which arc finally averaged to obtain the quantum fidelity, 
= ~M Sjli This averaging procedure al- 

lowed us to reduce considerably the level of statistical 
fluctuations. 

Results. — Our numerical results showed an unusual 
behavior for the quantum fidelity, which stays equal to 
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FIG. 1: Fidelity decay for Az = 3L ho , N A = 10 5 , and three 
values of the perturbation: e = 10 -3 , e = 10 -5 , e = 10~ 7 (the 
critical time increases with decreasing perturbation). The 
dashed lines are fits obtained using Eq. ([2} with ui z T = 4.86. 



unity until a critical time tc, and then drops rapidly to 
small values (Fig. 1). The critical time is defined as 
the time at which the fidelity has dropped to 60% of its 
maximum value, i.e. F(jc) — 0.6. Interestingly, the 
fidelity can be nicely fit by a Fermi-like curve 



/(*) = (1 - /oo) 



1 + exp 



t-r c 
T 



+ U 



(2) 



Equation ^ reveals the presence of two distinct time 
scales: (i) the critical time tq and (ii) T < Tc, which 
measures the rapidity of the fidelity decay. The param- 
eter /oo simply reflects the fact that the fidelity cannot 
decay to zero, because the system is confined in a finite 
region in space. In Fig. 1, u z T — 4.86 is the same for 
all three cases and is equal to the oscillation period of 
a particle trapped in the anharmonic potential V(z), for 
an initial condition z(0) = Az, i(0) = 0. This means 
that the fidelity decay occurs over one single oscillation 
period. 
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FIG. 2: Fidelity decay for iV A = 10 5 , e = 10" 5 , and different 
values of the displacement Az — 2Lho, Az — 3Lho, and Az = 
4Lho- The steeper curves correspond to larger values of Az. 
The dashed lines are fits obtained using Eq. ((2]), with uj z T = 
4.27 (Az = 4L ho ), lj z T = 4.86 (Az = 3L ho ), and lj z T = 5.50 
(Az = 2L ho ). 

This behavior was confirmed by investigating the de- 
pendence of the quantum fidelity on the initial displace- 
ment Az (Fig. 2). For each value of Az, the parameter 
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T appearing in Eq. ([2]) is taken to be equal to the corre- 
sponding oscillation period in the anharmonic potential 
V(z). The oscillation period decreases with increasing 
energy (and thus with increasing Az) and indeed the 
fidelity drop becomes steeper for larger displacements. 
The critical time also slightly increases with decreasing 
displacement and goes to infinity for Az — > 0. This pre- 
sumably happens because the dynamics of the conden- 
sate becomes too regular when the confinement is har- 
monic. 
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FIG. 3: Critical time tc (in units of uj^ 1 ) versus perturbation 
amplitude e, for Az = 3Lho and Na = 10 5 . The solid line 
represents the curve tc ~ — to Ine, with ui z to — 3.3. 

Figure 3 shows that tq depends logarithmically on 
the perturbation amplitude, i.e. tq ~ —to hi e, with 
uj z to ~ 3.3 (this is the straight line depicted in Fig. 3). 
This is similar to what was obtained for another self- 
consistent model [l9j], suggesting that such behavior is 
generic for TV-body systems, at least in the mean field 
approximation. The critical time also depends on the 
number of condensed atoms Na- By decreasing Na, tc 
becomes considerably longer (Fig. 4), and for Na — ► 
(i.e. for the linear Schrodinger equation) we have that 
Tc — » oo. For sufficiently large condensates (Na > 2xl0 4 
for the case of Fig. 4), the critical time depends logarith- 
mically on the number of atoms. 

Finally, the collapse of the quantum fidelity 
is clearly linked to the phases of the wavefunc- 
tions. Indeed, by defining an 'amplitude fidelity' 
Fa(t) = (J \ipH 1 ipH 2 \dx) (which neglects information 
on the phases), we have verified that F a (t) shows no sign 
of a sudden collapse when the ordinary fidelity drops. 

Discussion. — We have shown that the evolutions of 
two BECs in slightly different Hamiltonians diverge sud- 
denly after a critical time tq- The interaction between 
the atoms is obviously a vital ingredient, as tc — > oo 
when the coupling constant vanish, i.e. for the linear 
Schrodinger equation. The crucial point is that, for the 
GPE, the unperturbed Hamiltonian Hq depends on the 
wave function. When the perturbation induces a small 
change in -0, Hq is itself modified, which in turns affects 
tp, and so on. Thanks to such a nonlinear loop, the per- 
turbed and unperturbed solutions can diverge very fast. 
In contrast, for the single-particle dynamics Hq is fixed, 
and the solutions only diverge because of the perturba- 



250 
200 - 



150 - 




10 10 10 10 10 10 

N 4 

FIG. 4: Critical time tc (in units of uj^ 1 ) versus number of 
atoms Na in the condensate, for Az = 3Lho and e = 10~ 5 . 
The straight line is a guide to the eye. 



tion 6H. Changes in ip add incrementally to each other, 
but cannot trigger the nonlinear loop observed in the GP 
simulations [221 ]. 

This behavior is clearly linked to the phases of the 
wavefunctions and could be tested experimentally by 
studying the effect of a random potential on the inter- 
ference pattern of two condensates @, 0, H, Q ■ A possible 
experiment could be performed as follows (see Fig. 5). 
First, a BEC is created in a single- well trap; then, the 
trap is deformed into a double- well potential 9, 23], with 
the barrier between the wells sufficiently high that the 
two condensates cannot tunnel through it. The conden- 
sates are left in the double trap for a time long enough 
to reach their ground state. A laser speckle is then used 
to create a small random potential of amplitude e and 
correlation length a z . If o~ z <C d, where d is the distance 
between the two BECs, each condensate is subjected to a 
different random potential with the same statistical prop- 
erties. 

In order to excite the dynamics, the total double-well 
trap is suddenly shifted by a distance Az of the order of 
a few micrometers. The BECs evolve in their perturbed 
trap for a certain time t, after which both the trap and 
the random potential are switched off, so that the con- 
densates can overlap and interfere. We predict that the 
contrast of the interference fringes will depend on the 
time t and on the perturbation e, in a manner analogous 
to the quantum fidelity: if t < Tc(e), the contrast should 
be large, whereas it should drop significantly for times 
larger than tc ■ Performing several experiments with dif- 
ferent perturbation amplitudes and different numbers of 
condensed atoms should allow one to reproduce qualita- 
tively the logarithmic scalings of Figs. 3 and 4. Accurate 
time-resolved measurements might even reproduce the fi- 
delity drop time T. 

These results, together with those obtained for an elec- 
tron gas at low temperature [l9j . suggest that many- 
particle systems display a generic sudden decay of the 
quantum fidelity. Thanks to the ease with which ultra- 
cold atom gases can be created and manipulated in the 
laboratory, BECs should constitute an ideal arena to test 
these predictions experimentally. 
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FIG. 5: Two BECs (ipHi and V , h 2 ) are created in a double- well 
trap. A random potential with amplitude e and correlation 
length a z is superimposed to the confining potential. In order 
to excite the nonlinear dynamics, the trap is suddenly shifted 
by a distance Az. 
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